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Abstract 

The persistent current in three-dimensional mesoscopic rings is investigated 
numericahy. The model is tight-binding one with random site-energies and 
interaction between electrons. The self-consistent Hartree-Fock approxima- 
tion is adopted for the interaction between electrons, and models with up to 
6 X 6 X 20 sites are investigated. It is shown that the long-range Coulomb 
interaction enhances the persistent current for rings with finite width. The 
origin of the enhancement is discussed. 
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I. INTRODUCTION 



It has been predicted by Biittiker et al. that an equihbrium persistent current flows 
in a mesoscopic normal-metal ring threaded by a magnetic flux. This prediction has been 
experimentally confirmed 10-0. Especially, the experiment on single Au ring |Q revealed 
that the size of the persistent current observed was quite large. Namely the observed current 
was / = (0.3 ~ 2.0)/o, where Jq = evp/L is of the order of the persistent current in an ideal 
one-dimensional ring with vp being the Fermi velocity and L being the circumference of the 
ring. The current is periodic in the magnetic flux piercing the ring with period 0o = h/e 
as expected. Thus the persistent current was predicted theoretically and confirmed by the 
experiment. However, this is not the whole story. Quantitative understanding of the current 
has not been established. 

It has been clarified by numerical and theoretical investigations that theory 
without taking into account the interaction between electrons can give current which is too 
small compared to the experiments Inclusion of the mutual interaction into theory is 

hard and difficult to assess the validity of the approximation employed, but several groups 
PJ^ have done such calculations, and obtained results that sample averaged current is of the 
order of {le/L)!^, where le is the elastic mean free path. In this case the current is periodic 
in 0o/2, because sample averaging kills the component with period 0o- This current is much 
smaller than Jq in the experimental situation of L >> /g. In these theories the Coulomb 
interaction acts mainly to maintain the local charge neutrality of the system. Another 
theory by Kopietz [0] takes into account long-range Coulomb interaction energy associated 
with the charge fluctuation. This theory gives the average persistent current of the order of 
0.1 Jo with period 0o/2. This value is large enough to explain the experiments 0. However, 
since the fluctuation should be suppressed by the Coulomb interaction itself to maintain 
the approximate local charge neutrality, it is questionable that such large value survives 
improvement of the theory [[TI] , p!2[ . 

On the other hand, numerical investigation is free from approximations, and has a pos- 
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sibility to give a clear answer to the origin of the large persistent current. However, unfor- 
tunately rigorous treatments are limited to quite small system sizes. Thus until now only 
small one-dimensional systems have been investigated with negative answer to the role of 
the Coulomb interaction p!3|-p^. Namely, it has been clarified that in one-dimensional sys- 
tems the repulsive interaction between electrons suppress the persistent current in situations 
relevant to the experiment. We have also investigated the one-dimensional system with the 



Coulomb interaction between electrons . Our results were the same as others. However, 
we also found that the Hartree-Fock approximation gives almost the same results as those 
of the exact diagonalization. The merit of the Hartree-Fock approximation is that we can 
deal with much larger systems. Thus we now apply our method to three-dimensional rings 
to investigate the effect of the Coulomb interaction. In the present work we have found that 
for the three-dimensional ring the Coulomb interaction enhances the persistent current. 



II. MODEL 

We consider the Anderson model with Coulomb interaction between electrons. For sim- 
plicity we neglect spin freedom of electrons. Thus we consider a tight binding model on a 
cubic lattice. Our model ring consists of x x A^^ sites: circumference iV^a, width Ny^a, 
and height N^a, where a is the lattice constant. The coordinate of the i-th site is given by 
three integers {ii,Wi,hi), where I < ii < Ni, I < Wi < N^^, and 1 < hi < N^. Periodic 
boundary condition is imposed only in ii, namely £i + Ni = £{. In the other two directions 
the lattice is truncated. Each site has random site-energy Si, which has uniform distribution 
between —W/2 and W/2. Thus our Hamiltonian is 

id i 

+1 V—i4ci - Po)ic]cj - po), (1) 

where c|(q) is the creation (destruction) operator of a spinless electron, tij is non-zero and 
equal to t only for nearest neighbor hopping, 6ij is non-zero only for hopping along the 

3 



circumference of the ring and has a value 27r0/A^£0o, V = e'^/Anea is the strength of the 
Coulomb interaction, and po is neutralizing positive charge at each site. To make a ring 
which has the same lattice constant everywhere we have embedded our three-dimensional 
ring in the four- dimensional space, so the distance Vij between i, j sites is given by 

a i\£ 

+ {wi-Wjf + {hi-hjy, (2) 

where R = [2 sin(7r/A^f)]~"'^ is the radius of the ring. 



The following treatment is the same as before |16|: The interaction term is decoupled 
following standard Hartree-Fock scheme by introducing the average (cfc-,), which is deter- 
mined self-consistently by iteration. Once the ground state energy Eg{(f)) is determined 
self-consistently, the persistent current is given by the derivative of the energy Eg{(j)) 
with respect to the magnetic flux: 

Since Eg{(f)) is an even function of and periodic with period 0o, we first Fourier analyze 
the energy, 

Eg{4>) = Eo + Y.En cos(27rn-f ), (4) 

n=l ^0 

and obtain the Fourier components of the current: 

= E-i-^nSin(27rn^) 

n=l ^0 
°° t (J) 

= '^In—sin{27(n—). (5) 

n=l ^0 (po 

Here J„ is a dimensionless number, which depends on realization of the random site-energy 
and the strength of the Coulomb interaction. In the following we mainly concentrate on the 
00-periodic component, Ji. 
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III. RESULTS 



We have investigated up to 720-site systems by this method, the largest being Ni^ — 20, 
A^^„ = 6, and = 6 with 300 electrons. The randomness is chosen such that the sample is 
in the diffusive regime, which is realized around W/t — 4 for the samples we have considered. 

For the largest sized system we have investigated five samples, which mean five different 
realization of the random site energies. The results for /i for these five samples are shown in 
Fig. 1, where samples are distinguished by different symbols, circles, two kinds of triangles 
(normal and upside down), squares, and diamonds. The absolute values of Ii are plotted, 
the sign of Ii is distinguished by open and closed symbols, open meaning positive and closed 
meaning negative value of /i, respectively. It can be seen that the current is enhanced in 
four samples out of five. This is the first evidence that the Coulomb interaction can enhance 
the persistent current. 

On the other hand, it has been established that the Coulomb interaction suppresses the 
persistent current for one-dimensional samples. It is evident that the finite cross section of 
the sample is responsible for the enhancement. Thus to clarify the effect we have calculated 
persistent current for systems with the same length but with smaller cross sections. In such 
calculations we have chosen the number of electrons such that it is about 40% of the lattice 
sites. The results are shown in Fig. 2. Here we show sample averaged typical current, which 
is defined as /typ = ((^i))^^^, where ( ) means average over five realizations of the random 
site energies for each sample size. The system sizes considered are 1 x 1 x 20 sites with 8 
electrons, 2 x 2 x 20 sites with 32 electrons, 3 x 3 x 20 sites with 72 electrons, 4 x 4 x 20 
sites with 128 electrons, 5 x 5 x 20 sites with 200 electrons, and 6 x 6 x 20 sites with 300 
electrons, which are shown by open circles, closed circles, open triangles closed triangles, 
open squares, and closed squares, respectively. It is seen that the typical current in the 
absence of the interaction increases quite weakly as the cross section of the ring, namely the 
number of channels, increases. This is the feature expected in the diffusive regime. On the 
other hand, when the interaction is taken into account, the current begins to depend on the 
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cross section. The current is suppressed for the one-dimensional case, but is enhanced for 
the three dimensional case, the largest sized system being the strongest. 

In order to understand the origin of the enhancement, we have divided current into three 
parts, namely into the single-particle term Ig, the Hartree term 7h and the Fock term Ip, 
which are defined as follows, respectively: 

i^i i 

lH{<i>) = -id/d4>)J2Via/2r,,){clc,){c]c,), 

and 

lFi<P) = id/d4>)J2Via/2r,,){clc,){c]c,). 

The root mean square of these contributions for the one-dimensional system and the largest 
sample, 6 x 6 x 20 sites, are shown in Figs. 3(a) and (b), respectively. We see that in the 
one-dimensional case the contribution from the interaction terms are smaller than the single 
electron term. Is, which decreases as V increases. On the other hand both the single-particle 
term and the Hartree term are enhanced considerably by the Coulomb interaction for the 
three-dimensional system, the Fock term remaining quite small. The former two contribution 
has different signs, and the strong enhancement is almost canceled in the total current. 

IV. DISCUSSION 

In the present work we have shown for the first time that the Coulomb interaction 
can enhance the persistent current. It is true that the enhancement is not large enough 
to be comparable with the experiment, where the persistent current of the order of pure 
one-dimensional ring, Iq — evp/L, has been observed. The size of the current amounts 
to h — An/Ne ~ 0.628, if we use vp — 2\t\a?i, and N(, — 20. However, the cross section 
dependence of the enhancement is encouraging, since the actual sample has width of 90 nm, 
which is much larger than the lattice constant of the Au, 0.3 nm. 
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Then the next question will be whether the strength of the Coulomb interaction is appro- 
priate or not. This is a difficult question, since size of our system is still too small compared 
to the experiment, where about 10^ atoms are involved. Here we only comment that if our 
system is actually made of Au atoms, t is of the order of leV, and bare Coulomb energy 
between one lattice constant is larger than t: V — e^/Anea ~ 5eV > t. So even if we 
consider the optical dielectric constant e(oo) by ion core, V and t are of the same order. 
We have found the strongest enhancement when y is a httle smaller than t. In the actual 
system the interaction is screened by electrons, but such screening effect is included in our 
self- consistent solution. We need the flux dependence of the total ground state energy which 
includes the energy associated with the screening. Therefore use of the long-range Coulomb 
interaction is needed and justified. 

Our treatment is not exact as we use the Hartree-Fock approximation. We have justified 
it by comparison with the results of exact diagonalization for small one-dimensional systems. 
Actually for one-dimensional system the suppression of the current is stronger in the Hartree- 
Fock approximation than the exact method. Thus it is possible that the enhancement in 
the 3-d case is stronger in the exact treatment. However, more investigation will be needed 
in order for us to truly justify and assess our method. 

Now the final question will be the mechanism of the enhancement. One possible explana- 
tion is the screening of the random potential. The electron at the Fermi level feels screened 
site-energy which is fiattened, so the current is enhanced. In the case of one-dimensional 
system the screening is also effective, but electrons are not allowed to pass through each 
other, thus the current is suppressed by localization of some of electrons. When the cross 
section of the ring is small enough, this blockade is still effective. However, Figs. 3(a) and 
(b) suggest another scenario. It is seen that Ig and Jh are enhanced by the Coulomb interac- 
tion considerably for the tree-dimensional system. Similar behavior is also seen for systems 
with cross section larger or equal to 3 x 3 sites, but not seen for systems with smaller cross 
sections. The large value of the Hartree contribution means that the direct Coulomb energy 
changes considerably as the magnetic flux changes. This is caused by the deformation of the 
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electron wave function as the flux changes. The deformation of the wave function does not 
be done freely. The Coulomb energy suppresses it for the system to approach to the charge 
neutrality. Thus both the single particle part and the Hartree part are enhanced, but they 
nearly cancel each other. What left of this cancellation brings the enhanced total current. 
The finite size of the cross section makes the Coulomb energy larger, and also gives more 
freedom for the wave function deform as the flux changes. This makes the qualitative differ- 
ence between the one-dimensional system. Fig. 3(a), and three-dimensional system. Fig. 3(b). 
Figure 3(b) shows that the ground state is quite different from that of the non- interacting 
system, where we get quite small persistent current. This explains why the theories using 
the diagrammatic expansion PJ^JT^|T^ could not get the enhancement: Since the ground 



state has changed considerably from that without the Coulomb interaction, perturbational 
calculation cannot give correct answer. 

In conclusion we have shown for the first time that the persistent current is enhanced 
by the long-range Coulomb interaction for three-dimensional rings. We have clarified that 
a ring with finite cross section is necessary for the enhancement. 
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FIGURES 

FIG. 1. Fourier component of the persistent current with period is plotted as a function 

of the strength of the Coulomb interaction V jt for the randomness W jt = 4.0. The five different 
symbols, circle, two kinds of triangles, square, and diamond indicate five different realization of 
the randomness. The absolute values of the current are plotted with the open (closed) symbols 
indicating that the value is positive (negative). 

FIG. 2. Root-mean-square of the Fourier component of the persistent current with period (f)o, 
ityp, is plotted as a function of the strength of the Coulomb interaction V/t for the randomness 
W/t = 4. Results of six different sized systems are plotted. The system sizes considered are 1 x 1 x 20 
sites with 8 electrons, 2 x 2 x 20 sites with 32 electrons, 3 x 3 x 20 sites with 72 electrons, 4 x 4 x 20 
sites with 128 electrons, 5 x 5 x 20 sites with 200 electrons, and 6 x 6 x 20 sites with 300 electrons, 
which are shown by open circles, closed circles, open triangles, closed triangles, open squares, and 
closed squares, respectively. 

FIG. 3. Root-mean-square of the Fourier component of the three contributions to the total 
current for the systems with (a) 1 x 1 x 20 sites with 8 electrons and (b) 6 x 6 x 20 sites with 
300 electrons. The single-particle contribution Jg, the Hartree contribution In, and the Fock 
contribution /p are shown by circles, triangles and squares, respectively. It should be noticed that 
the vertical scale is expanded 10 times in (a). 
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